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Abstract. We develop a continuum theory of linear viscoelastic response in oriented monodomain nematic 
elastomers. The expression for dissipation function is analogous to the Leslie-Ericksen version of anisotropic 
nematic viscosity; we propose the relations between the anisotropic rubber moduli and new viscous coeffi- 
cients. A new dimensionless number is introduced, which describes the relative magnitude of viscous and 
rubber-elastic torques. In an elastic medium with an independently mobile internal degree of freedom, the 
nematic director with its own relaxation dynamics, the model shows a dramatic decrease in the dynamic 
modulus in certain deformation geometries. The degree to which the storage modulus does not altogether 
drop to zero is shown to be both dependent on frequency and to be proportional to the semi-softness, 
the non-ideality of a nematic network. We consider the most interesting geometry for the implementation 
of the theory, calculating the dynamic response to an imposed simple shear and making predictions for 
effective moduli and (exceptionally high) loss factors. 
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1 Introduction 



The equilibrium mechanical response of liquid crystalline 
elastomers can be soft or hard depending on the relation 
between the imposed strains and the nematic director, 
in particular, if the director is able to respond by rotat- 
ing. With such unusual equilibrium elasticity one might 
expect dynamical response to be equally unusual. If the 
elastic forces are small, then the return to equilibrium is 
driven more weakly than in conventional systems. How 
does the dynamics of internal director rotation, and the 
corresponding time-dependent softening of rubber-elastic 
response, determine the dynamic mechanical response of 
a nematic rubber to a small amplitude oscillatory shear? 
In a dynamic-mechanical study of Gallani et al. G[ a 
' stress response to an imposed oscillating shear has been 
studied. Although a polydomain elastomer has been ex- 
amined, the authors reached a conclusion that the re- 
sponse "is insensitive to the isotropic-nematic transition" 
and only obtained a non-trivial result in the smectic-A 
phase. Subsequent studies of aligned monodomain elas- 
tomers [|j also did not find any unusual effect in the ne- 
matic phase and went on to investigate the mechanical 
effects in the smectic-A phase. The elastic properties of 
smectic rubbers are very interesting on their own, with a 
number of spectacular effects even in equilibrium, see |3|,^], 
^,^|. However, here we would like to address a more basic 
and physically clear case of nematic elastomers, that is, 
rubbery networks with an aligned uniaxial anisotropy of 
their polymer strands. 



Some of the reasons that no exceptional effects were 
found in the nematic phase by |j],|2| could be that (i) the 
nematic region for the materials studied was only ~ 7° 
and smectic pre-transitional effects were important, (ii) 
the authors aimed to plot the whole range of dynamic 
modulus, including very high glassy values at low tem- 
peratures, thus masking a subtle nematic region and (iii) 
they did not study low enough frequencies. More recently 
[Q, it has been demonstrated, both theoretically and ex- 
perimentally, that a dramatic reduction of storage modu- 
lus G and the associated increase in the loss factor tan S 
should be expected and indeed found in monodomain ne- 
matic elastomers sheared in certain geometries (cf. Fig. ^|). 
This effect allows one to directly probe the basic equilib- 
rium properties of nematic rubbers and also access the 
new kinetic parameters - viscous coefficients and relax- 
ation times. 

In order to study the dynamics of mechanical response, 
one needs to model viscoelastic properties, that is, de- 
scribe viscous dissipation in a system moving towards its 
equilibrium. At present, there is no microscopic model 
that would even approximately describe the dynamics of 
anisotropic rubbery networks. However, as in the contin- 
uum theory of liquid crystals, much progress can be made 
on a phenomenological level, using the symmetries of vari- 
ables contributing to the physical effects. 

The continuum dynamics of anisotropic fluids describes 
macroscopic phenomena on the level of coarse-graining 
analogous to that of the nematic Frank elastic free energy. 
By movements of a "liquid particle" one understands bod- 
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ily translations and also changes of local director orienta- 
tion of a physically infinitesimal volume, generally of the 
correlation size £n ~ lOnm, including many molecules in 
thermal equilibrium with the reservoir, and characterised 
by the local nematic tensor order parameter Qij(T) = 
Q(n;nj — 5<5jj). Two physical fields describe the state of 
motion of nematic liquid: the local variation of director 
orientation 6n(r,t) = n — n with respect to the equi- 
librium n and the fluid velocity v(r,t), the time deriva- 
tive of the corresponding elastic displacement u(r,t) in 
the description of elastic response. The original deriva- 
tion of hydrodynamic equations for the nematic liquid, 
due to Ericksen and Leslie is presented in some detail in 
key monographs on liquid crystals, e.g. [^J. Here we only 
discuss the matters relevant for the basic description of 
nematic elastomers, omitting many fine and subtle points 
of this complicated subject, e.g. ||. We develop a for- 
malism describing the linear viscoelasticity of an elastic 
medium with an independently relaxing director degree of 
freedom and make a connection between the (anisotropic) 
rubber-elastic and viscous coefficients in such a system. 
The nematic dynamics coupled to the underlying rubber 
elasticity leads to an expression for the entropy production 
density, the Raylcigh dissipation function in the Leslie- 
Ericksen form, determined by the fluxes: i - the symmet- 
ric strain rate and 4?{f2 — u)) - the rate of relative rotation 
between the nematic director and the elastic matrix: 

Ts = Ai{n ■ g ■ nf + 2A 4 [n x g x nf + 4A 5 ([n x g ■ n]f 
+ I 7l [nx|(fi- W )] 2 + 72n.i.[nx|(fi- W )] , 

where the new linear viscous coefficients, A{, are linear 
combinations of the classical Leslie coefficients, which in 
elastomers can take extraordinarily high values com- 
pared with simple nematics. 71 and 72 describe the same 
kind of losses involving the rotating director as they do in a 
conventional nematic liquid. We find a direct proportional- 
ity between the symmetry-related elastic and viscous con- 
stants, connected by the corresponding relaxation times 
for each deformation mode. One of the main results of 
our analysis is a new dimensionless number, which we call 
Ne = rjv/(L^,), where r\ is a viscosity, /i is a rubber shear 
modulus and L is a characteristic length. Ne determines 
the local balance between viscous and elastic torques in 
the material, in analogy with the Ericksen number Er || . 

An essential part of description developed in this paper 
is the separation of time scales. We argue, on the basis of 
comparison with available experiments and qualitative ar- 
guments, that the characteristic relaxation time of director 
motion is much greater than the time scale of mechanical 
relaxation in an ideal polymer network. Accordingly, we 
only consider the low-frequency limit of rubber viscoelas- 
tic response, reflected in essentially linear frequency de- 
pendence of complex modulus G* = G + iujr/. In contrast, 
the much slower director relaxation allows us to examine 
both its low- and high-frequency regimes. 

This paper is organised as following. The next Sec- 
tion, continuing the introduction into the subject, briefly 
sketches the fundamentals of equilibrium nematic rub- 



ber elasticity, focusing on the linear continuum descrip- 
tion rather than the full molecular theory valid to high 
deformations. Following this, in Section [3| we revisit the 
Leslie-Ericksen formulation of the continuum dissipation 
function Ts and show that, as would be expected by sym- 
metry, it is identical to the elastic energy density only 
with strains substituted by strain rates. We then discuss 
the fundamentals of nematic rubber viscoelasticity. In Sec- 
tion 4 the practical calculation of dynamic moduli G* (to) is 
presented for three principal shear geometries, as an exam- 
ple of general theory implementation. The results indicate 
which combinations of elastic and viscous coefficients are 
relevant for the response and, therefore, directly measur- 
able by experiment. We then conclude by discussing the 
role of soft elasticity in slowing the mechanical relaxation 
and make contact with new experimental data. 

2 Nematic rubber elasticity 

Equilibrium elastic properties of monodomain nematic rub- 
bers are well-studied, both theoretically and experimen- 
tall y a nd are described at some length in review articles 
[[n)],[ll[] . A full molecular theory of ideal nematic networks 
gives the elastic free energy density 

F=±c x k B TTr(± T -£ e 1 .±-Q + §B(Det[AJ - l) 2 , (I) 

where A is a Cauchy strain tensor, A;k = dRjdR^. In 
a nematic elastomer network, £ are the uniaxial matri- 
ces of chain step-lengths before (0) and after the director 
n has rotated by a certain angle 9 during the deforma- 
tion: £ij = £j_Sij + [£|| — l^]nirij. The last term in ([!]), 
the bulk-modulus contribution independent of the con- 
figurational entropy of polymer chains, is determined by 
molecular forces resisting the compression of a molecu- 
lar liquid, B ~ 10 9 J/m 3 , much greater than the typical 
value of rubber modulus fi ~ cJt^T, with c x the effective 
crosslinking density, giving fi ~ 10 5 J/m' 3 . This large bulk 
energy penalty constrains the value of the strain determi- 
nant, Det[A] ~ 1 (which, in other words, means that the 
material is physically incompressible when subjected to 
all strains except hydrostatic compressions and dilations). 
Eq. (|l|) is only limited by the Gaussian network assump- 
tion and is valid up to very large strains where mechanical 
softness, rotational instabilities and optical effects occur - 
well beyond the limits of continuum theory. 

Olmsted Jl^] , in studying soft elasticity - shape changes 
in nematic elastomers without energy cost, has derived 
the small-deformation continuum limit of the full expres- 
sion (|]), when Ay = 5% + diUj and |Vu| <C 1. We shall 
use a slightly different notation, more suitable for the 
linear-elastic description because it better complies with 
the standard textbook formalism of uniaxial elasticity, e.g. 
Jl3| . The small-deformation limit of ([[]) is 

F = Ci{n-g-nf + 2C 2 Tr[|](n • ~g ■ n) + C 3 (Tr[g]) 2 
+ 2C* 4 [n x £ x nf +4C* 5 ([n x|-n]) 2 (2) 
+ \Di[n x(!2- a;)] 2 + D 2 n ■ g ■ [n x {fi - u)]. 
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where n is the undistorted nematic director, = — 
•jTrfe] Sik is the traceless part of linear symmetric strain 
£ik = ^(dkUi + diUk), which is the only variable of classical 
continuum elasticity |k| . In a system with an internal ori- 
entational degree of freedom, the nematic director with its 
own relaxation and dynamics, the antisymmetric part of 
strain expressed by the local rotation vector fl = icurlw 
may contribute to the physical properties. Analogously, 
the small vector u> is a convenient measure of director ro- 
tations, u) — [n x 5n]. In fact, it is the relative rotation, 
the difference nx(f2 — a>), that causes the elastic response 
and a number of effects unique to nematic elastomers |14f . 

One expects that in a rubber or dense polymer melt 
the bulk modulus C3 is very large, C3 ~ B. We, there- 
fore, shall only consider deformations with no bulk com- 
pression: Tr[s] = 0. In general, all other constants in the 
expression (||) are of the same order of magnitude, similar 
to the rubber modulus fx. The molecular model of ideal 
nematic elastomer gives specific forms for these con- 
stants: 

c x = 2C 4 = c K k B T, c 5 = | Cx fc B r (r ^ 1)2 , 

0-1)2 1 _ r 2 

D 1 = c x k B T± '-, D 2 = c x k B T , (3) 

r r 

in which case the condition for ideal soft elasticity holds, 



c 5 R 



c 5 



3 



= 0. 



(4) 



Model expressions for elastic constants (|J) depend, apart 
from the universal rubber-elastic energy scale [i, on a sin- 
gle parameter r. In the molecular model ([!]) of ideal ne- 
matic polymer networks one finds that r = i\\/£±, the 
ratio of the principal step lengths of the anisotropic poly- 
mer backbone (or equivalently r = (R\\/R±) 2 in terms of 
the principal values of gyration radii tensor) . In non-ideal 
elastomers, this parameter is more complex, determined 
by a number of other factors, for instance when there are 
fluctuations in composition, see the Appendix. Neverthe- 
less, it has to be a function of nematic order parameter 
Q(T), satisfying a linear limit r w 1+/3Q, at least at small 
Q. In the isotropic phase, at Q — and r = 1, the elastic 
constants (||) become, as expected: C\ — 2C4 = 2C5 = /i, 
Di = D 2 = and the elastic energy (^) reduces to a stan- 
dard Lame expression. See the Appendix for a discussion 
of non- ideality, that is where Crf ^ 0. 

In an incompressible material, all deformations are es- 
sentially shears, albeit sometimes seen in a rotated frame. 
Figure [l] shows the character of these shears in the eq. (||) 
and their corresponding elastic constants. C\ "lives" along 
the director, C4 in the plane perpendicular to no in which 
properties are isotropic (ekk, £mm and e m k are of the same 
status since the solid is uniaxial). Most interesting are C5, 
involving shears such as e n k, which span directions parallel 
and perpendicular to the initial director and which can in- 
duce it to rotate; Di which resists director rotations with 
respect to the solid matrix; and D 2 where the rotation is 
coupled to the symmetric shears also involved in C5. The 




Fig. 1. The symmetries of the shears and rotations possible in 
a uniaxial solid at constant volume and with an internal ori- 
entational degree of freedom. The modes are labeled with the 
appropriate elastic constant and symmetric shear or rotation. 
C4 also shows a section perpendicular to the director to illus- 
trate shears e m k not involving n (the directions m and k are 
perpendicular to n). In the two coupling terms, Di and D2, 
the small director rotation « |<5n|. 



inter-relation between these three processes is what gives 
rise to the effect of soft elasticity - here shape change of 
the type e n k without energy cost. If the director is allowed 
to react ideally and is regarded as a slave variable, then the 
rcnormalised resistance to deformations e n k, that is C§, 
vanishes as we see in eq. (^). We shall see, in Section 4 
for simple shear geometries, the mechanism by which this 
is achieved - at least at harmonic order (since these sim- 
ple geometries do not allow the full shear and extensional 
freedom required for softness at all amplitudes). 



3 Nematohydrodynamics of elastic solids 

When nematic elastomers are strained, there is director ro- 
tation along with stresses - the latter yielding both body 
forces and torques. At finite strain and rotational rates 
there will be stresses of both elastic and viscous origin, the 
latter being conceptually parallel to those arising in clas- 
sical liquid nematics. We shall accordingly review classical 
nematohydrodynamics, ignoring the subtleties that reside 
at the heart of the subject. There are important differ- 
ences with liquids - if torques are induced by flow, they 
will be balanced in solids, on length scales longer than the 
nematic penetration depth, by torques generated by the 
matrix, that is, by Di and D 2 elastic terms that involve 
the anti-symmetric part of the strain A. The issue of possi- 
ble director gradients and corresponding Frank elasticity 
in nematic elastomers has been extensively discussed in 
the literature It is known that, unless there 

are special reasons for a director singularity (such as in 
disclinations or narrow domain walls), Frank elastic effects 
play a minor role in the free energy balance and can be 
neglected. Instead, the rubber-elastic matrix generates the 
torques, which shape the final character of the dynamics. 
A new dimensionless number will be introduced to replace 
the Erickscn number describing a torque balance in ordi- 
nary nematics. There is another important difference - in 
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liquids, as flow proceeds, the strain and the director ro- 
tation can increase without limit, whereas in solids both 
are constrained (and are coupled to each other) by the 
matrix. Elastomers are capable of huge distortions but we 
shall limit ourselves to linear theory. Nematohydrodynam- 
ics in the elastically highly non-linear regime is a subject 
of another study. 



3.1 Leslie-Ericksen formulation of nematic viscosity 

The equation of motion can be written in the usual form 
of local balance of forces: in vector components, 



at 



(5) 



where p is the density and ilki the tensor of momentum 
flux density, 77^ = pv^Vi — a"ki+-P<i>ki with P the local pres- 
sure and o"ki the stress tensor. In an isotropic Newtonian 
liquid <7ki is proportional to the symmetric strain rate, 
^4ki = | (c^k^i + <9i«k) = £ki- In uniaxial nematic liquid 
crystals, the stress tensor depends not only on the fluid ve- 
locity gradients, but also on the components and the gradi- 
ents of the local nematic order Qki ■ Because rotations and 
corresponding torques are involved if the director changes 
differently from the local fluid rotation, the stress tensor is 
no longer symmetric and also depends on the relative ro- 
tation combinations N — -^8n+ i [n x curl v] || , which is 
nothing but a time derivative of the relative-rotation com- 
bination, that is N = [n x ^(17 — a;)], cf. equation (||). 
As a result, Leslie writes the viscous stress tensor 

cry = airiirij^kmft-k^m + a 2 n- 1 N- i + a^Ni (6) 
+aiAij + a^riin^A^j + aen^n^A^ 

(assuming incompressible fluid, Ay^ = divv = 0). Its anti- 
symmetric part contributes to the local torque r = [nxh] 
(cf. ||) where h is the molecular field, h derives also in 
part from Frank and external field terms. These contribu- 
tions are balanced by the anti-symmetric terms deriving 
from viscous flow processes: 



h 9 ( 

">i = "o 



SF 



9F 
orii 



+ 72«k-4ki (7) 



with 7i = ay, — a.2 and 72 = 0:3 + a-i = ae — a§ (the 
Parodi relation, a representation of the Onsager principle 
of kinetic coefficient symmetry). Here the viscous coeffi- 
cients a±, ... , ae all depend on the magnitude of nematic 
order parameter Q(T). One can show that in an ordinary 
nematic liquid crystal near the weak first order transition, 
as Q —* 0, they should behave as |lq| 



ai cx Q 2 , (a 2 , as, a 5 , a 6 ) cx Q 
also 71 cx Q 2 and 72 cx Q. 



(8) 



Thus, in the isotropic phase, only one of the Leslie coef- 
ficients survives: «4 — > 2r). Molecular theory (l7|[l|| also 
shows that in a typical nematic liquid of rod-like molecules 



these coefficients may have very different magnitude and 
even sign: a%, 03, ae - negative, 04, as - positive; (|aa| ~ 
05) ^> |cK3 1 , \ae\; a\ is generally small and may be pos- 
itive or negative in different materials; far from the ne- 
matic transition (at Q — > 1) the "isotropic coefficient" 
«4 ~ |ct2 1 ~ a§. De Gennes assembles the experimental 
values for Leslie coefficients in MBBA at 25°C ||, while 
the monograph by de Jeu [p^[ provides the values for an- 
other classical nematic, PAA at 122°C: in 10 -3 Pa.s 



MBBA 
PAA 



ai a 2 «3 
6.5 -77.5 -1.2 
4 -6.9 -0.2 



a 4 a 5 ae 71 72 
83.2 46.3 -34.4 76.3 -78.7 
6.8 5 -2 6.8 -7.1 



3.2 Nematic elastomers 

We now develop Leslie-Erickson theory for solid nematics 
by a direct analogy. The nematic dynamics is coupled to 
the underlying anisotropic elasticity described in its linear 
limit by eq. (0) and sketched in Fig. |l|. It is thus natural to 
develop the dynamics using the same symmetry-grouping 
of terms as in (||), rather than as in (|^) as is usual for 
classical nematics. 

Two differential equations, (||) with the viscous stress 
given by eqs. (JsJ) , and (Q), form the complete set describing 
the linear viscous effects in the nematic fluid. Neglecting 
the effects of heat convection, the total energy dissipation 
(the entropy production) in such anisotropic medium is 
expressed by the volume integral of the conjugate forces 
and fluxes: 



TS = J dV (oijAj + (h ■ n) - (F • !?)) 

= J dVfaiij + ih-N)). (9) 

We re-write the density of dissipation function in a form 
matching the elastic energy density (Q): 

Ts = Ai(n-g- n) 2 + 2A 4 [n x g x n] 2 + 4A 5 ([n x | • n]) 2 
+ ±7i-/V 2 + 72" ■£■ N (10) 

in the fully incompressible case. Thus, by differentiation of 
( |l0| ) respectively by the symmetric strain rate e\j and by 
iVj, one obtains a representation of the symmetric viscous 
stress tensor and the nematic molecular field, contributing 
to the local torque, analogous to the parallel expressions 
(||) and @ for simple nematics: 

cry = 2Ai (n ■ e • n) ruiij + 4A4 [n x (n x i x n) x n] . . 



+ 4A 5 ^ [ (n x e • n) xnj.rij 

+ I72 {niNj + Nitij) ; 
hi = 71 + 72«j£ij , 



[( 



n x e • n I x n . n. 



(11) 



where ofj is the symmetric part of an expression like (|6j) 
which emerges naturally since we express Ts in terms of 
symmetric and anti-symmetric variables separately. Here 
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the constants are linear combinations of classical Leslie 
coefficients 



A 5 = \{2a 4 + a 5 + a 6 ). 



(12) 



In the isotropic limit, one finds A\ = 2A 4 = 2A§ — > 77, 
reminiscent of the isotropic Lame limit of the correspond- 
ing elastic constants. 



3.3 Relaxation times 

The direct correspondence between the elastic constants of 
the energy density F, eqs. (|j), and the viscous coefficients 
of the dissipation function Ts, eq. (|l0|), gives: 



Ai -> d; 71 -> L>i and 72 



Do 



Of course, this is not a surprise since the symmetry of 
an equilibrium elastic deformation and of a viscous flow 
in a uniaxial continuum is the same. Equally, there is a 
direct correspondence in the dependence of the various 
coefficients on nematic order parameter Q, see eq. (||) for 
viscosities and eq. (||) for elastic constants. We can equiva- 
lently, therefore, represent the viscous coefficients of a ne- 
matic elastomer by products of the corresponding rubber- 
elastic constant and an appropriate relaxation time, that 
is, 

Ai = CiT R ; 71 = Dm; 72 = D 2 t 2 . (13) 

In principle, the relaxation times for the various Aj will 
be different, and different from the Tj. However, one ex- 
pects a significant separation between the two groups of 
characteristic time scales, that is between the trs and the 
TjS. The director rotation time, t\ — Ji/Di, has been 
used in a simplified, single- viscosity analysis and then 
estimated experimentally as T\ ~ 10 _1 — 10~ 2 s. In con- 
trast, the characteristic time of rubber relaxation tr is 
expected to be much shorter. A recent statistical theory 
explicitly calculated this time from the analysis of net- 
work strand relaxation within a tube model p0( ] , showing 
this is of order of Rouse time. One could also confirm this 
by a following qualitative observation: The linear viscous 
stress (11 ) and the analogous corresponding expression for 
the elastic stress could be combined, producing the effec- 
tive frequency-dependent moduli in a form (C; + Sb)A\). 
This is the low-frequency limit of a general complex mod- 
ulus G*(uj), showing the rubber plateau modulus G = C\ 
and the initial rise in the loss modulus with frequency, 
describing a viscous flow. In the classical picture of poly- 
mer dynamics, e.g. pM , the next characteristic regime is 
at a Rouse frequency, when the signal cannot propagate 
along the polymer chain length and mechanical response 
is provided by individual segments, thus causing a climb 
of G*(oj) towards the glass plateau. Therefore, the char- 
acteristic time scale tr in the estimation relations ([l3]) is 
of the order of the Rouse time, tr ~ 10~ 4 — 10~ 6 s. 

Could there be viscous softness analogously to the clas- 
tic softness of ideal elastomers, for instance in eq. (Q)? Two 



related arguments suggest that this is not possible. The 
viscous combination A§ — 7! /871 , analogous to the renor- 
malised elastic constant , can be rewritten with the 
benefit of timescales eq. ( |l3|) as: 



7I/871 



mn 8£>i 



(14) 



Even if we took the rotational relaxation times equal, 
Ti = T2, the translational relaxation time tr would ap- 
pear to upset the possibility of renormalising (jlj) to zero, 
even if we had — > 0. Elastic softness arises because an 
anisotropic distribution of chains can be rotated undis- 
torted and thus at constant entropy. However, individ- 
ual chains will be distorted and there must reasonably 
be dissipation associated with their flow relative to the 
matrix. Thus elastically soft distortions should have asso- 
ciated dissipation, except perhaps accidentally if rates of 
71 , 72 and A$ relaxations compensate to make (O) vanish. 



3.4 Symmetries and order parameter 

The ability to neglect complicated incrtial effects in fluid 
dynamics is controlled by the small magnitude of Reynolds 
number, Re = pvL/rj <c 1, with L the characteristic 
length and 77 the typical viscosity. In nematic liquids, an- 
other dimensional number is introduced to characterise 
the relative magnitude of hydrodynamic and Frank elas- 
tic torques contributing to eq. (g), where the left hand 
side is ~ K\7 2 n, while the right hand side is ~ jVv, 
with the rotational viscosity 7 ~ r/. Small Ericksen num- 
ber Er = rjv L / K -C 1 means that the director orienta- 
tion is mostly controlled by equilibrium elastic free energy, 
while at Er ^> 1 the director generally follows the lo- 
cal orientation provided by the flow. In elastomers we are 
concerned with the balance between flow-induced torques, 
again scaling as ~ ijVii, and those of the rubbery matrix, 
expressed by dF/d9 ~ D\Q. This yields a new dimension- 
less group characterised by the number Ne = r/v/(LDi). 
Note that the assumed domination of rubber-elastic over 
Frank effects essentially means that D\ 3> K/L 2 , that 
is, Ne -c Er. Our analysis of an example of simple shear 
deformation, in the next Section, spans the full range of 
small to large Ne, which we shall define more precisely for 
that geometry. 

We have seen in Section 2, in the discussion of eqs. (||) 
and chain-anisotropy parameter r, how the anisotropic 
rubber moduli depend on the underlying nematic order 
parameter Q, usually a function of temperature or sol- 
vent concentration. In the isotropic phase, at Q = and 
r = 1, the elastic constants return to the classical Lame 
values Ci = 2C 4 = 2C* 5 = fi, £>i = D 2 = 0. The cou- 
pling constant D 2 depends on the linear power of Q. As 
a result, when there is anisotropy, Q ^ and r ^ 1, 
the sign of the elastic constant D 2 depends on whether 
r > 1 (prolate order, D 2 < 0) or r < 1 (oblate order, 
D 2 > 0) . The sign of director rotation relative to the ma- 
trix, uj — fl, varies accordingly to whether chains are pro- 
late (they align with the extension direction associated 
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with the shear, see Fig. [l] for £> 2 ) or oblate (alignment 
along the compression axis). In any event, a term such as 
-D2£nkW, bilinear in e n k^, always reduces the elastic en- 
ergy, irrespective of the sign of D 2 , and indeed this is the 
actual mechanism whereby a shape change can ideally be 
achieved with no energy cost jl2],[2^] . There is an appealing 
analogy between this behaviour and the Leslie coefficient 
72 of discotic nematics. It is known (24|,^1 that the dif- 
ferent (disk- instead of rod-like) molecular shape leads to 
the opposite sign of 72, with the according consequences 
for the flow alignment properties. 

Ideal rubber-elastic expressions (||) suggest that the 
coupling constant D\ is proportional to (r — l) 2 ~ Q 2 . 
Continuing the analogy with liquid nematic viscosity, we 
should recall the estimates of Imura and Okano [|6| , but 
also a following qualitative consideration. The rotational 
viscosity 71 is determined by the antisymmetric part of 
viscous stress tensor Oy , which cannot be proportional to 
the linear power of the symmetric tensor nematic order 
parameter Qy, but at least its square. However, there is 
a delicate problem arising in the analysis of soft elastic- 
ity. The renormalised shear modulus C§ = C5 — D 2 /8Di 
should reduce to the bare C5 — ► ^/j, in the isotropic phase, 
at Q — ► 0. However, the ideal values of D\ and D2 re- 
sult in a finite renormalisation in this limit. One must 
revise the conclusions of non-ideal, semi-soft theory |p"lf . 
Eq. ( |30|) and the discussion in the Appendix suggest that, 
in fact, Di has a additional small (semi-soft) correction 
~ Q. This resolves the problem of making the renormal- 
isation D2/8D1 vanish at Q — > 0, but raises a question 
about the symmetry consideration that 71 (and D±, by 
parallel) cannot be linear functions of Qy. The paradox 
is safely resolved when one recognises that the semi-soft 
coefficient ai in cq. ( |30| ) is a linear function of Qf, the or- 
der at network formation. Therefore, in fact, both D\ and 
71 in nematic elastomers depend on the bilinear combina- 
tion Qjj^Qij, which means ~ \Q\, and thus no symmetry 
problem arises. In an ordinary liquid nematic there is no 
issue of formation order being frozen by crosslinking, and 
7i ~ Q 2 as expected. 



3.5 Balance of forces and torques 

To complete the general analysis of this Section, we briefly 
discuss the stresses and torques that follow from the com- 
bination of elastic and dissipation functions. In an over- 
damped system, traditionally ignoring inertial effects at 
low-frequency, we have two equations of motion: the bal- 
ance of forces and torques. The first condition requires 
locally balancing the total symmetric stress tensor 



OF 
~de 



d(Ts) 
di 



(15) 



In the case when the relaxation time scales are separated, 
as in the discussion above, the viscous stress contribution 
to the force balance is minor. In contrast, the balance of 
torques, r = [n x h] = in the absence of external fields, 



requires the full molecular field 

h = h e \ + h visc = D\ [n x (fi — u?)] + D 2 n ■ e + (16) 
+ 71 [n X (fl - w)] + 7 2 n • g . 

Here both groups of terms are manifestly of the same order 
of magnitude in the regime of frequencies characterising 
director rotations. We shall see in the particular calcula- 
tion of Section 4 that the condition of zero local torque 
allows one to obtain, for instance, the rate of director vari- 
ation 9. 

In the limit of isotropic rubber Q — > the only relevant 
equation is that for the stress, which in this case reduces 
to 

a_ = C±e + A^i , 

which is a simplest viscoelastic approximation at lowest 
frequencies, or shear rates, of a general linear-response ex- 
pression a(t) = J G(t — t')e(t')dt' . Another known limiting 
case is that of an uncrosslinked nematic liquid crystal. In 
this case all equilibrium (zero-frequency) elastic moduli 
Ci and Di are zero and the only contribution to the stress 
and torque balance are the Leslie-Ericksen eqs. (|TT[). In a 
low molecular weight nematic, this is the full description 
of an anisotropic Newtonian liquid. In a polymer nematic, 
again, one expects a complex viscoelastic response func- 
tion, G(t — t'), with several characteristic time scales, from 
the shortest Rouse time, to the entanglement and diffu- 
sion times (if applicable) |2l| . Equations (|ll|) are thus the 
low-frequency limit of such complex nematic viscoelastic- 
ity. 

The theory of elastomer nematohydrodynamics has 5 
phenomenological viscous coefficients. In a simplified single- 
viscosity model we have found remarkable effects the 
internal director relaxation has on the macroscopic dynamic- 
mechanical response of a nematic rubber. Here we follow 
the example of ordinary liquid crystals, where both the 
simplified and the full Leslie-Ericksen formalism have been 
used successfully over the years. In the practical calcula- 
tion implementing the above ideas, we now move from the 
one-constant analysis of M to a full analysis. 



4 Shear stress and its relaxation 
4.1 Simple shear deformations 

In a study of linear response, we shall first examine three 
principal simple shear geometries, as shown in Fig. |^. 
These are also the geometries that one achieves in a typical 
dynamic- mechanical experiment . A geometry of uniax- 
ial extension, more commonly found in studies of equilib- 
rium stress-strain in elastomers, is less appropriate for an 
oscillating regime because of possible slow relaxation [^6| 
and incomplete sample recovery on each cycle. The sim- 
ple shear e(i), assumed externally applied to the sample 
(Fig. ||), is the single xz-component of the full Cauchy 
strain, the same for each director setup G, D and V. It 
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Fig. 2. The geometry of simple shear experiment with three 
principal orientations of the nematic director no, labelled G 
(for no along the shear gradient), D (displacement) and V 
(vorticity). The small-amplitude shear e M ~ \ee}^ 1 is applied 
to the elastomer and the measured stress <j(uj) provides the 
linear response modulus in each of the three configurations. 



also automatically satisfies the necessary incompressibil- 

(s) 1 

ity constraint. We have then the symmetric part e\x = 2 £ 
and the antisymmetric (body rotation) part expressed by 

Note that the three principal shear geometries in Fig. g 
are the same as in the classical setting for Miesovicz viscos- 
ity experiments. That is, if the director is kept immobile 
(e.g. by a strong external field, which incidentally would be 
hardly possible in elastomeric network under strain) , then 
the orientations V, D and G correspond to the Mieso- 
vicz viscosities r] a , rj^ and rj c , respectively. Of course, the 
whole point of the present paper is to examine the effect 
the director freedom to rotate is having on the viscoelastic 
response. 

We have argued above that Frank effects are subordi- 
nate in the overall balance of torques. Perhaps, this ap- 
proximation needs to be reconsidered if a sheared sample 
is very thin in the z-direction: if the director is anchored 
on top and bottom surfaces (which are the plates in a 
shear-sandwich experiment), one expects a Frank energy 
density of the order K/d 2 . The elastic energy density is 
set by ix and the two scales take on equal importance when 
d = £ = \J Kj [i, the nematic penetration depth. Taking 
typical values K ~ 10 _11 J/m 3 and fi ~ 10 5 J/m 3 one has 
£ ~ 10~ 8 m. Thus in practical situations with a sheared 
sample thickness d ~ 100 Lira, and more, ignoring the ef- 
fect nematic director gradients seems to be a safe first 
approximation. 

The remaining free energy density, eq. (^]), takes the 
form, in the three cases of Fig. ||: 



= (C 5 + 

F B = (C B + 
Fy = C4 e 2 



±[D 1 -2D 2 ])e 2 -±{D 1 -D 2 )e6 
£[£>i + 2D 2 }) e 2 - |(Z>i + D 2 )e 9 - 



(17) 



where the small change in director orientation, Sn, is taken 
equal to the angle 9. Clearly, one does not expect director 
rotation to occur in the "log-rolling" geometry V. 

The bilinear term ~ e 9 clearly reduces the energy: if 
a strain e is imposed then the director responds by ad- 
justing 9 to lower energy. Given sufficient time to reach 
equilibrium, 9 adopts its optimal value for a given defor- 
mation e. Returning this minimised value $g,d in (18) to 



@, the free energy at a given strain is also optimal 
Di T D 2 



>G,T> 



2D! 



Ft 



G,D 



c 5 



n 2 , 
2 * e 2 (18) 



8D 



Notice that if the nematic elastomer is ideal, then the free 
energies in the cases G and D vanish (C5 — D 2 /8Di = 0) 
because their geometry allows the director to respond to 
the shear and internally relax. Case V remains elastically 
hard. In fact, because of the chosen restricted strain ge- 
ometry, Fig. ||, the response of even ideal elastomers is ac- 
tually quartic, rather than completely soft (true softness 
requires some unconstrained extension as well as shear). 
The molecular model fll],^2| yields the quartic penalty 
F = 2 A 1 (r-1) 2 > wn i c h we are neglecting in the present 
linear-response analysis. 

The generalised force driving the director rotation an- 
gle 9 if) to equilibrium is in fact a torque, dF/d9. In F it is 
the externally imposed shear e(t), with its time variation, 
that is the ultimate driving agent. The flatness (softness) 
of the free energy F(e) will make the resulting dynamical 
response characteristically slow in the geometries where 
director rotation is possible E^j. Alternatively, if the im- 
posed dynamics, e.g. e(t), is fast compared with that of 
0(t), then one will not attain the ideal states (18) and the 
moduli will become frequency-dependent and deviate, for 
dynamic reasons, away from softness. It is this rheological 
subtlety that we now wish to examine. 



4.2 Viscoelastic response 

Driving a nematic elastomer by imposing a component of 
strain, in general, leaves open the possibility of dynami- 
cal response of not only the director but also the other 
components of strain (and thus involving for instance the 
translational viscosities A\, A$ and A5). Then there will 
be several coupled dynamical equations for the system, 
see for instance ^7J where the response of a nematic elas- 
tomer to a step-extensional strain is calculated. Here we 
confine ourselves to the dynamical response to imposed 
simple shear with other components of strain clamped. 
An elastic torque acting on nematic director n in cer- 
tain deformation geometries is resisted by a linear viscous 
torque. Describing the director rotation by a small angle 
6, cf. Fig. ||, and continuing to ignore inertial effects, the 
balance of torques and forces is expressed by 



5F_ 

59 



SiTs) 



= 0. 



(19) 



The dynamical equation describing the evolution of the 
director approaching its equilibrium is given, for the two 
geometries where director rotation is present: 

(G :) 11 6 = -D l e + ±{D 1 -D 2 )e{t) 

+|(7i- 72) i{t) (20) 

(D :) j 1 e = -D 1 9 + l(D 1 + D 2 )e(t) 
+ |(7i+ 72) e(t) • 
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These linear, inhomogeneous differential equations are eas- 
ily solved. After the transient relaxation 9 ~ 9^e"^ Dl ^ 11 ^ t 
associated with starting the strain oscillations has com- 
pletely relaxed, the steady-state response is given by the 
particular solutions with e ~ e luJt : 



0g,d(*) - / 

J — O 



dt > e -(D l/71 [t-t')) 



(21) 



x I ^eM + ^T 72/7i)^') 



where the signs — or + correspond to the G or D geome- 
try, respectively. The solutions depend on a characteristic 
time for director relaxation, t\ = 71/-D1, or equivalently, 
the characteristic frequency of the response scales like 
-D1/71. We can now give a more concrete expression for 
the new dimensionless number Ne = r yiVv/(Di9), since 
9 ~ e and Vw ~ e ~ we, whence: 



Ne = jiU>/Di = Tito 



(22) 



Below we shall see that all depends on t\U) - when it is 
small, we have soft or semi-soft equilibrium elastic re- 
sponse, when it is large, we have hardening because di- 
rector relaxation does not keep pace with changing strain. 

The linear "nominal stress" in response to the imposed 
simple shear deformation e(t) is given by the sum of elastic 
and viscous stress functions er = cr e i ' 



<r' visc , that is 



dF 



d(Ts) 
de 



(23) 



In the principal shear geometries, the stress expressions at 
a given frequency of imposed strain take the form 

ctgM = 2 (C 5 + \[D X - 2D 2 }) e(u) - |(A - D 2 ) 9{w) 
+2(A 5 + i[ 7l - 2 72 ]) e_i( 7l _ 72 )0 (24) 

a D (cj) = 2(C 5 + |[A + 2D 2 ]) e(w) - |(A + £> 2 ) #H 
+ 2(A 5 + i[ 7l + 2 72 ]) £-i( 7 i+72)0 (25) 

ctv(w) = 2C 4 e(w) + 2^ 4 e • (26) 

One finds two kinds of viscous stress terms, of different 
orders of magnitude. The contribution 7$ ~ 7 e is of the 
same order as the elastic terms, as indicated by the torque 
balance eq. (|20|). In contrast, the terms ~ A^i are of the 
order Ae ~ C{lut^)e <C Ce at frequencies below Rouse 
values. The linear viscoelastic theory we are considering 
is applicable at much lower frequencies, where the most 
interesting physics is due to director relaxation and me- 
chanical softness. 

The effective response modulus in the "log-rolling" ge- 
ometry V is unchanged by the nematic director dynam- 
ics, G(ui) = 2C4 + 2iwAi (the loss modulus negligible, as 
discussed above). In two other geometries, where the di- 
rector rotation does take place, the dynamic modulus is 
modified by the internal director relaxation. Substituting 
the Fourier transforms of equations ( |2l"| ) into the expres- 
sions (p4| ) and (p5|), we obtain the nominal stress in the 
form o(lu) = GJlu)e(lu). Remarkably, although perhaps 





Reduced frequency cot 



Storage modulus G'/u 



(b) 



Fig. 3. (a) The frequency dependence of storage (G' , circles) 
and loss (G" , squares) moduli in the two shear geometries G 
and D in units of rubber modulus /i. There is no loss in the V- 
geometry and, since G' v = /i, the plot represents the variation 
of the ratios G' G /G' V and G' D /G' V . At u -> both G' -» 2Cf . 
(b) The Cole-Cole plot, a variation of G" with G' . A semicir- 
cular shape indicates a single relaxation time response. In both 
plots, solid lines are for Q = 0.58, dashed lines for Q = 0.35. 



predictably, the response in these two geometries is exactly 
the same: ctq = od, despite the difference in the rotations 
9q and (9 D . The corresponding storage and loss moduli 
are given by the real and imaginary parts of the effective 
complex modulus G(u>). It has a single-relaxation time be- 
haviour with a characteristic frequency oj\ = D\ / 7 i or the 
corresponding relaxation time t\ = 1/uj\. The dynamic 
moduli, in both G and D geometries, can then be written 
in a universal form: 



G'(u) = 2(C 5 - Dl/8Dx) 



(loti) (Am - -D172) 



(-D271 ~ -D172) 2 



. + (^n) 2 4D l7l 2 
un (^i(72/7i) 2 ~2A 5 /n) 



(27) 



(28) 



Figure ^| shows the example of frequency dependence 
of eqs. ( p7| ) and (|2^). In order to plot these functions, we 
need to make several numerical assumptions about the 
parameters. We thus take, arbitrarily for illustration pur- 
poses, r — 1 + 1.5Q at a fixed temperature with values of 
the nematic order parameter corresponding to deep in a 
nematic phase, Q — 0.58 (so that the chains are weakly 
prolate, r w 1.87), or corresponding to near a transition 
point, Q — 0.35. The ratio 72/71 ~ 1 and the viscous 
constant A$ ks O.IT1C5. Finally, the semi-soft addition to 
the coupling constant D\ (see the Appendix) is taken as 
a\ = 0.1. For these values, plots of G and G" reveal the 
expected single-relaxation behaviour. At zero frequency of 
imposed strain oscillations the real (storage) modulus in 
both shear geometries G and D is equal to 



G 



1 _ o/^R 



2C 5 H = 2(C 5 - D5/8D1) 



The measure of non- vanishing renormalised shear modu- 
lus Cf- is the characteristic parameter of semi-softness 
in non-ideal nematic elastomers. The high frequency re- 
sponse for lot\ ^> 1 is 



G'oo = 2C 5 - i A (72/71) + |^i(72/7i) 2 



(29) 



Terentjev and Warner: Viscoelasticity of nematic elastomers 



9 



The characteristic single-relaxation time behaviour seen 
in fig. H is due to our assumption about separation of time 
scales, T\ 3> tr. For that reason, we are able to distinguish 
the plateau in G , the maximum and the "high-frequency" 
drop in G" and a semicircular shape of Cole-Cole plot 
p{| . For the chosen value of A$ , one can begin to see the 
next rise in G" , when the frequency approaches the next 
characteristic point, the Rouse frequency leading towards 
the glass plateau value G' ~ 10 11 Pa at high frequency. 
That transition would correspond to another semi-circular 
Cole- Cole plot at much higher values of G'; we can see 
the beginning of that graph segment in fig. ||(b). With the 
assumed separation of nematic and rubber relaxation time 
scales (and the assumed ideal polymer network, with no 
entanglements and free dangling ends, which are known 
to lead to slow relaxation), the mechanical losses at low 
frequencies are only through the lag in director rotation. 

The initial expressions, eqs. ( p4| ) and (]25|), differ in 
sign of terms with D2 and 72 constants. This reflects the 
tendency to align or repel the director by these terms, 
the roles being reversed when we interchange prolate and 
oblate symmetry (r > 1 and r < 1, the sign of both D2 
and 72 reversing on this exchange). However, when the 
director relaxation (|2l]) is substituted back into the stress 
expressions, one only finds the products -D272, and a corre- 
sponding unique effective complex modulus. Another fea- 
ture of G*(lo) is the non-dimensional ratio 72/71, which 
is a familiar parameter in the dynamics of ordinary liquid 
crystals || in the context of director flow alignment or 
tumbling. As the brief discussion of modified torque bal- 
ance in Section 3 indicates, one does not expect a steady- 
state tumbling in elastomers, where there is memory of the 
original director through the elastic energy. Tumbling and 
analogous effects may only be found in transient regimes 
or in non-linear elasticity, to which we return elsewhere. 

Variation of dynamic complex moduli with tempera- 
ture is very interesting to examine graphically, as well as 
analytically. Figure ^(a) shows the dependence of storage 
modulus G' on the reduced temperature T/T„i, for a se- 
quence of increasing frequency values. In order to produce 
these plots, we must take a further assumption about the 
variation of the order parameter Q(T). We take, rather 
arbitrarily, Q » (1 — T/T ni ) ' 33 simply because this was 
an approximate fit to the experimental measurement of 
Q(T) in M. We also take, following classical nematics, 

71 = giQ &n<i 72 = 32Q and further fix 32/31 = \ to have 

72 ~ 7! at Q r~j 0.5. These rather limiting assumptions 
nevertheless usefully illustrate the qualitative behaviour 
of the linear elastic response functions both deep in the 
nematic phase and near the assumed "critical point" T nl . 

The characteristic time of director relaxation t\ = 
71/-D1 is a function of temperature through its dependence 
on nematic order parameter. Both the coupling constant 
D\ and the rotational viscous coefficient 71 are functions 
of Q(T), as discussed in Section 3. Hence the relaxation 
time n should have a weak Q-dependence in the nematic 
phase, which should change to t± ~ \Q\ in the vicinity of 
the nematic- isotropic transition point. We discuss subtle 
limit problems and the form of G(T, u>) as a function of T 
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Reduced temperature T/T 
(a) 



Reduced temperature T/T 



(b) 



Fig. 4. (a) Plot of reduced storage moduli G'/fi against re- 
duced temperature for a number of increasing reduced frequen- 
cies: Co = 0.1 (circles), 0.5 (squares), 1 (triangles), 2 (diamonds) 
and 10 (stars). Moduli approach the semi-soft constant G' a at 
zero frequency and low temperature, cf. Fig. [](a). An apparent 
critical behaviour near the transition point T n j is due to the ne- 
matic order parameter variation Q(T) appearing in the elastic 
and viscous coefficients, (b) Plot of loss factor tan 8 = G" /G' 
against reduced temperature for a number of increasing fre- 
quencies. The curves are labeled the same as in (a). 



for various fixed u) in the Appendix. The combination u)T\ 
can be rewritten as 



71 
1 



Q 2 



fij oi|Q| + (r-l) 2 /r' 



where we recall that r ~ 1 + 0G (in practice, the linear 
relation holds to high values of Q |7],[31j). We shall continue 
taking — 1.5 in our illustration. One can also scale w 
by n/gi to give a dimensionless frequency ui; thus wt\ = 
Luf(Q) where the non-dimensional function f(Q) is clearly 
/~|Q|asQ^0. 

The cusps near T n j seen in Fig. ||(a) are a consequence 
of our assumption of the critical form Q(T) ~ (1— T/T ni ) 33 . 
For a real, non-ideal monodomain nematic elastomer with 
of necessity a high-temperature paranematic state and a 
residual small Q above T n ;, the cusps are rounded off but 
the general form survives, including the non-monotonic 
variation G'{T) at high frequencies. We note in passing, 
that the experimental result for G'(T) is very close to 
these predictions: although the very high frequency was 
not achieved, the gradual departure of G'(T)/fx from the 
universal low-frequency curve was clearly registered. 

Figure |](b) also shows the unusually high loss tangent 
arising from the real and imaginary parts of complex mod- 
ulus. The new loss mechanism, due to internal director 
relaxation, is additional to all classical losses occurring in 
polymer systems - in particular those occurring near the 
glass transition. The information about the latter is, as be- 
fore, contained in the essential frequency and temperature 
dependence of rubber-elastic constant fi, of which the ex- 
pressions for F and Ts are the low- frequency limit. There 
is a striking contrast between the fig. 0(b) with unusually 
high values of tan S over a broad temperature range and 
the traditional distinct loss peaks in classical polymers 
and rubbers. Physically, the maximal loss occurs where 
the the imposed strain frequency, w, matches the director 
relaxation rate 1/ri . We have a crossover from the direc- 
tor keeping up with strain and therefore allowing the the 
system to be very soft, to where the director fails to re- 
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spond quickly enough and the system becomes hard. The 
large change of modulus when one crosses this frequency 
interval is manifested, by the Kramers-Kronig relations, 
in a very high loss peak. 



5 Conclusions 

We have developed a model of time-dependent response 
and rheology of nematic elastomers by coupling nematic 
rubber elasticity with nematohydrodynamics. The theory 
represents the low-frequency limit of general viscoelastic 
description of nematic polymer networks (which is not 
known). By making use of significant separation between 
the natural nematic relaxation time scale and that of a 
rubber network, we examine the dissipation at relatively 
low frequencies due to internal director relaxation in a 
deformed elastomer. Because of the lagging of director ro- 
tation in response to time-dependent imposed shears, soft 
or semi-soft response can only be partially achieved, the 
hardening depending on frequency and on temperature. 
The nemato-elastic-hydrodynamic equations are derived 
in generality and solved for the rheology of a particular, 
important example, that of simple oscillating shear, giv- 
ing rise to unusual temperature dependence of response 
moduli and loss tangents (tan 5 ~ 1.5 — 2). Thus these re- 
sults will be important for a wide variety of applications in 
damping, acoustics, frequency and directional dependent 
materials, and so on. 

We have already reported on rheological experiments 
|0] that show the qualitative features of our theory. These 
and light scattering experiments |2j],|3(| give values of 
Frank constants, rubber constants C5, D\, D% and the 
rotational viscosity coefficient 71 for nematic elastomers 
(and the 71 also for the corresponding polymer melts). 
The C's and D's can also be deduced from several static 
measurements, in all cases giving typical rubber values 
~ 10 3 — 10 5 Pa, modulated by the discussed Q-dependence. 
Frank constants take values typical of classical nematic 
liquids, K ~ 10 -11 N. But the 7 values are enhanced by 
polymerisation. This is not surprising since on reorienting 
a chain one must translate it bodily. An extended object 
sweeps through space and the dynamics are complex. In a 
network this is even less straightforward because the chain 
is mechanically tied to a matrix and must resolve its mo- 
tions with those of other chains. It turns out that the 7's 
are further enhanced, reaching values of order ~ 10 3 Pa.s 
and more, leading to the mentioned separation of charac- 
teristic time scales. 

We appreciate valuable discussions with S.M Clarke, H. 
Finkclmann, T.C. Lubensky, R.B. Meyer, P.D. Olmsted 
and W. Stille. 



A Semi- softness and rheology 

At fixed temperature in the nematic state, Fig. || shows 
the rheological consequences of nematohydrodynamics in 



elastomers - as frequency drops and director relaxation 
had keep pace with imposed strains, the response gets ever 
softer. In the limit ujt\ — > it either vanishes (ideal soft- 
ness) or becomes very small (semi-softness) as has been 
confirmed by many quasi-static experiments. We wish to 
comment on the mechanisms for softness and non-ideality 
and then how, at various fixed frequencies, one can see 
how G tends to conventional rubber elastic response as T 
increases and Q — » 0. 

Soft shear response occurs because anisotropic network 
chains accommodate a macroscopic shape change by ro- 
tating their distribution of shapes at constant entropy. 
When chains tend to isotropy, Q — ► and r — ► 1, shape 
change can only be accommodated by chains distorting 
their distribution, thereby decreasing their entropy and 
increasing their free energy. The modulus is then /i. We 
must see how this classical limit is achieved; after all the 
cancellation C$ = C5 — D 2 /SDi = appears to hold for 
all r including r — * 1. 

Semi-softness is expected when there are fluctuations 
of composition, rod-like crosslinks, sources of random ne- 
matic field, and any other form of non-ideality that pro- 
hibits the finding of an isotropic reference state p^p3|. 
Then there are additions to the Trace formula (Q) which 
are of the form a sin 2 8 for simple shears, see for instance 
[p2| . The degree of semi-softness, a, can be calculated di- 
rectly from various models of non-ideality, or measured. 
For continuum theory such additions appear in the eqs. (|^) 
for the Ci and Di, and upset the cancellations in (||), yield- 
ing a Ck* ^ 0, which is the non zero limit at u)t\ — » 
in Fig. for G'. However, the additions inspired by a 
clearly must vanish as Q — ► and r — > 1 and the deli- 
cate limit question must still be resolved. In particular in 



(| Cf - C 5 



j/i as r ■ 



In the ideal case there is no high temperature order, 
softness is perfect and for T < T n \ the modulus C;f is 
identically zero. Several systems, depending on their ther- 
momechanical history, are close to this [[22|. Since ideally 
Q jumps to zero at T = T n j, there is no limit problem. D\ 
and D 2 cease to exist, the analysis ([l8]) is invalid and C5 
is not renormalised. 

In the semi-soft case, (Q) places bounds on the form 
the non ideal additions to D\ can take. Since D\ ~ Q 2 
and D 2 ~ Q, in the ideal case the limit D 2 /Dx is finite 
as Q — > 0. If however Di has additions D\ ~ (Q 2 + a\Q), 
then the limit of D 2 /D\ vanishes as Q — > 0, eliminat- 
ing any renormalisation. Explicit calculations agree with 
symmetry arguments that the additions to D± are indeed 
~ Q and reflect the thermomechanical history of the mate- 
rial. For instance fluctuations (due to compositional fluc- 
tuations in the polymers that make up the network) in 
the effective order felt by a chain, Q e g = (1 + S)Q (with 
(5) = 0), yield for the modulus: 



Dx = n' 



{r-lf 



ail 



(30) 
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where // 



dete 



dcte 



/j, and 



_q / ^ Q/(3 + 2Q/) 

[ l (i + 2Q f )(i-Q f y 



Here formation conditions for the network are denoted by 
a subscript /: the step length tensor at formation is £ 
with order parameter Q / . The determinant-factors reflect 
spontaneous shape changes since formation which can be 
very large (up to several hundred %). The sign of the coef- 
ficient a\ is that of Q f . It is the order at network formation 
that induces the residual order Q at high temperatures 
and the sign of that order follows that of the formation 
order. Thus the combination a±Q is always positive and 
thus the non-ideal additions to D\ in (|3C|) are really of the 
form | Q | , which we have used in section p.4| . 

Thus at sufficiently low frequencies one observes a rather 
small modulus G' = Crf > 0, which rises with increasing 
temperature to reach the classical value fi at Q — > at 
T n i. If, at fixed frequency on dropping the temperature 
the combination lot\ becomes appreciable (due to increas- 
ing n), the modulus G'(T) will depart from the universal 
low-frequency form G^ Tl ^ , see Fig. 0. 
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